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1. Introduction 

Traveling waves in nonlinear Hamiltonian lattice systems are ubiquitous in many 
branches of sciences and their mathematical analysis has attracted a lot of interest 
over the last two decades. In the simplest case of a spatially one-dimensional lattice 
with nearest-neighbor interactions - often called Fermi-Pasta-Ulam or FPU-type chain 
- the analytical problem consists of finding a positive wave-speed parameter a along 
with a distance profile R and a velocity profile V such that 

R'{x) — V(x + |) — V(x - |) , 


( 1 ) 
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is satisfied for all i6l. Here $ is the nonlinear interaction potential and the position 
Uj(t) of particle j can be recovered by 

Uj (t) = U(j-y/at), U(x) f V(y)dy, 

J xo 

which implies the identities 

iij(t) = \fcrV( K j — \/at) and Wj+i(t) — Uj{t) — R(j + | — \fat) 

for the atomic velocities and distances, respectively. In particular, u satisfies Newton’s 
law of motion 

Uj(t) = &(u j+ 1 (t) - Uj{t )) - &{uj(t) - Uj-i(t )), jeZ. (2) 

The existence of several types of traveling wave solutions (with periodic, solitary, 
front-like, or even more complex profile functions) can be established in different 
frameworks; see, for instance, [5, 2] for constrained optimization problems, [9] for critical 
point techniques, [8] for spatial dynamics, and [10, 13] for almost explicit solutions. 
However, very little is known about the uniqueness of the solutions to the advance- 
delay differential equation (1) or their dynamical stability within (2). The only nonlinear 
exceptions are the Toda chain - which is completely integrable, see [11] and references 
therein - and the Korteweg-de Vries (KdV) limit of solitary waves in chains with so called 
hardening. The latter has been investigated by Friesecke and Pego in a series of four 
seminal papers starting with [4], In this limit, solitary waves have small amplitudes, 
carry low energy, and are spread over a huge number of lattice sites. The discrete 
difference operators in (1) can therefore be approximated by continuous differential 
operators and the asymptotic properties are governed by the KdV equation, which is 
completely integrable and well understood. 

Another interesting asymptotic regime concerns solitary waves with high energy in 
chains with rapidly increasing potential. Here the profile functions localize completely 
since V converges - maybe after some affine rescaling - to the indicator function of 
an interval, see [3, 12] or [6] for potentials $ that posses a singularity or grow super- 
polynomially, respectively. The physical interpretation of the high-energy limit is that 
the particles interact asymptotically as in the hard-sphere limit, that means by elastic 
collisions only. 

The high-energy limit is another natural candidate for tackling the analytical 
problems concerning the uniqueness and the stability of traveling wave. In this context 
we are especially interested in the spectral properties of the linearized traveling waves 
equation - see the discussion in section 3.6 - but the convergence results from the 
aforementioned papers do not give any control in this direction. They are too weak and 
provide neither an explicit leading order formula for a nor the next-to-leading order 
corrections to the asymptotic profile functions. In this paper we derive such formulas 
and present a refined asymptotic analysis of the high-energy limit for potentials with 
sufficiently strong singularity. 
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In order to keep the presentation as simple as possible, we restrict our considerations to 
the example potential 

x 1 ( 1 

<P r) =-- -- rwr — mr — 

mn(m + 1 ) \ (1 — r) 

which is convex and well-defined for r < 1, satisfies 

<E>(0) = <h'(0) = 0, <h"(0) = 1, 

and becomes singular as r /* 1. The condition m > 1 is quite essential and shows up 
several times in our proofs. The other details are less important and our asymptotic 
approach can hence be generalized to the case 

is convex and smooth on some interval [a, b] with <E>'(a) = 0 
such that the limit lim^. ^ $(x )(6 — x) m does exist. 

This class also includes - after a reflection with respect to the distance variable - all 
Lennard-Jones-type potential, which blow up on the left of the global minimum. 

To simplify the exposition further, we merely postulate the existence of a family 
of solitary waves with certain properties but sketch in section 1.4 how our assumption 
can be justified rigorously. Specifically, we rely on the following standing assumption, 
where unimodal profile means increasing and decreasing for negative and positive x, 
respectively. 

Assumption 1 (family of high-energy waves). (Vs, Rs, cr <s)o<< 5 <i a f am ^y of solitary 
waves with the following properties: 

(i) Vs and Rs belong to L 2 (M) D BC^M) and are nonnegative, even, and unimodal. 

(ii) Vs is normalized by |||| 2 = 1 — 5 and Rs takes values in [0, 1). 

Moreover, the potential energy explodes in the sense of ps := f R $(Rs(x)) dx —> +00 as 
5 —>• 0 . 

Beside of S there exist two other small quantities, namely 

£5 ■— 1 — ^< 5 ( 0 ), ps := \J<Ts£™ +2 , (4) 

which feature prominently in the asymptotic analysis. The amplitude parameter £5 
quantifies the impact of the singularity and appears naturally in many of the estimates 
derived below. The parameter ps, which looks rather artificial at a first glance, is also 
very important as it determines the length scale for the leading order corrections to the 
asymptotic profile functions Vo and i?o- 

For the interaction potential (3) and the waves from Assumption 1, the existing 
results for the limit <5 —>- 0 are illustrated in Figures 1, 2 and can be summarized as 
follows. 


with m G ( 1 , 00) 


( 3 ) 
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Figure 1. Numerical examples of solitary waves for the potential (3) and as in 
Assumption 1: The graphs of Vg (black, dashed) and Rg (gray, solid) are plotted 
for m = 1.5 (top row) and m = 2.5 (bottom row). In the high-energy limit S —> 0 
(from left to right column), Vg and Rg approach the indicator function Vq and the tent 
map R 0 , respectively. See section 1.4 for details concerning the numerical scheme. 




Figure 2. Parameter plots for three different choices of m and the simulations from 
Figure 1: pg represents the potential energy, eg is the amplitude parameter from (4), 
and ag := ag/ (jj 2 £™) measures the relative deviation of the speed parameter a with 
respect to the asymptotic value from (32) 


Theorem 2 (localization theorem). In the high-energy limit, we have 

11^5 — Vo || 2 + ||-R<5 — -Rolloo + £ 5 + Rs - 1 0 

with 

Vq{x) := x( x ) Ro{ x ) '■= max{0, 1 — |x| } , 

where x denotes the indicator function of the interval [—+|], 

We give a short proof in section 1.3. The corresponding results in [3, 12] also 
provide lower and upper bounds but no explicit expansion for a§. 

1.2. Statement of the asymptotic result 

Our strategy for deriving a refined asymptotic analysis is to blow up the profile 
functions near the critical spatial positions and to identify equations that determine 
the asymptotic wave shape with respect to a rescaled space variable x. Specifically, 
we use the transition scaling in order to describe the asymptotic velocity profile near 
x — ±|, while the distance profile can be rescaled at both the tip position x = 0 and the 
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Figure 3. Schematic representation of the different scalings: The transition scaling 
describes the jump-like behavior of Vg near x = ±4 while the foot and the tip scaling 
magnify the turns of Rg at x ~ ±1 and x ss 0, respectively. The width parameter yg 
is introduced in (4) and satisfies yg ~ eg according to Corollary 13. 




Figure 4. Graph of the function So and its second derivative for m = TOi (black) and 
m = m 2 (dark gray) and m = m 3 (light gray) with mi < m 2 < m 3 . 


foot positions x = ±1, see Figure 3 for an illustration. Our main findings can informally 
be summarized as follows. 

Main results. In the high-energy limit h —)• 0, all relevant information on (Rs, Vs, as) 
can be obtained from the function So, which is defined by the ODE initial value problem 

§XW = yzr' 7—. s„(o) = s' 0 ( 0) = 0 (5) 

m + 1 (1 + S0U)) 

and plotted in Figure f. More precisely, 

(i) the velocity profile Vs converges under the transition scaling, 

(ii) the distance profile Rs converges under both the tip scaling and the foot scaling, 
(in) the rescaled parameters S m as, h -1 ^, and S^ps converge, 

where the respective limit objects can be expressed in terms of Sq and all error terms are 
at most of order 0(5 m ). 

The details concerning the convergence under the tip, the transition, and the foot 
scaling are presented in the Theorems 7, 9, and 15, respectively, while Corollary 13 
provides the explicit scaling laws for £s, ps, and as- Moreover, the combination of 
all partial estimates gives rise to the global approximation results in Theorem 16 and 
Corollary 17. 

The above results provide an improved understanding of the high-energy limit 
of solitary waves. In particular, it seems that our asymptotic formulas can be used 
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to control the spectrum of the linearized traveling wave equation, see the brief and 
preliminary discussion in section 3.6. 

The paper is organized as follows. In the remainder of the introduction we prove 
the localization theorem and discuss both the results from [3, 12] and the justification 
of Assumption 1 in greater detail. Then section 2 is devoted to the tip scaling, which 
turns out to be most fundamental step in our asymptotic analysis. In particular, we 
identify the intrinsic scaling parameters in section 2.1 and link afterwards in section 2.2 
the rescaled distance profile to the initial value problem (5). In section 3 we finally 
employ the results on the tip scaling and establish all other asymptotic formulas. 


1.3. Preliminaries 


In this section we prove Theorem 2 since it provides the starting point for our asymptotic 
analysis in section 2 and section 3. To this end it is convenient to reformulate the 
advance-delay-differential equation (1) as 

R = AV , aV = A&(R), (6) 

where the operator A stands for the convolution with the indicator function y. This 
reads 

i 


(AV)(x) = 


rx+7 


V{y) dy 


and the elimination of R reveals that (1) can be viewed as a symmetric but nonlinear and 
nonlocal eigenvalue problem for the eigenvalue a and the eigenfunction V. The proof 
that (6) implies (1) is straight forward and involves only differentiation with respect to 
x; for the reversed statement one has to eliminate the constants of integration by the 
decay condition V G L 2 (M). 

Using elementary analysis such as Holder’s inequality we readily verify the estimates 

||AU|| 2 < ||U|| 2 , HAUIU < ||U|| 2 , \\(AV)'\\ 2 < 2 ||U|| 2 , (7) 

for any function V G L J (M), and this implies that the potential energy 

V(V) := [ <h((AU)(x))dx, (8) 

Jr 

is well defined as long as ||U|| 2 < 1. Moreover, we get 

\\Rs\\oo = Rsi 0) = l-eg< IIVilla = 1-5 (9) 

for the family from Assumption 1. 

Lemma 3 (variant of the localization theorem). The estimates 

\\Vs-x\\ 2 <Ce s , \\Rg - AxWoo < Ce s (10) 

and 

c£ S < fig < Cy/T 6 (11) 

hold for some constants c, C independent of 5. Moreover, we have Es —> 0 as S —» 0. 
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Proof. We start with the identities 

RsiP) = {V s , X ), \\Vs-x\\l = m\l + 


2 {V s , x ), 


( 12 ) 


where (•, •) denotes the usual inner product in L 2 (M), and observe that (6) implies 

a s {l-S) 2 = (^(Rs),Rs). (13) 

Since (9) yields 5 < eg, we find 

o < II V s - xlll = (1 - 5f + 1 - 2(1 - E S ) < Ceg 

and hence (10)i, which in turn implies ( 10)2 thanks to Rg — A\ = A(Vg — x) and ( 7 ) 2 . 
By (6) we also have 


PjV«(0)=e; 


m+2 


$W(a;)) < ~ £s) < Ce 5 , 


and this provides the upper bound in (11) since the unimodality of Vg combined with 
(10)i guarantees that liminf^o V£(0) > 0. To obtain the corresponding lower bound, 
we notice that (1) ensures 

4||<E>W|U , C 


I R 


.//1 

■a 1100 


< 


(?8 


— _ _m+l 

crg£ s 


and hence 


_ _m+2 _ ,, 

&s £ g — ^5 


Rg(x) > 1 — Ceg for all |x| < 
due to Rg(0) = 0 and (0) = 1 — Eg. Combining this with (13) we obtain 


j4_ 

„ra+2 


(l-5) 2 > 




&(Rg(x))Rg(x) dx > c 


P5 


' -Us 


£ 


m +1 5 


and the proof of (11) is complete. Finally, the properties of $ imply 
Ps = V(V s )<£j m \\Rg\\l<C£j m 
so £g — y 0 is a consequence of pg —> 00 . 


□ 


1-4■ Justification of Assumption 1 

We briefly sketch how Assumption 1 can be justified using a constrained optimization 
approach. All key arguments are presented in [6] for non-singular potentials $ but can 
easily be adapted to the potential (3). At the end of this section we also discuss the 
results from [12] and [3]. 

The variational approach from [6] is based on the potential energy functional (8), 
which is convex and Gateaux differentiable on the open unit ball in L 2 (M); the derivative 
of V is given by 

d v V(V) = A&(AV ), 
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so the traveling wave equation (6) is equivalent to aV = dyV(V). We further introduce 
the cone C of all L 2 -functions that are even, unimodal and nonnegative, i.e. we set 

C := {V E L 2 (M) : 0 < V(x) < V(y) = V(—y) for almost all x < y < 0} . 


The key observation is that solitary waves as in Assumption 1 can be constructed as 
solutions to the constrained optimization problem 


Maximize V under the norm constraint ||V ^||2 = 1—5 and 
the shape constraint V E C. 


(14) 


In the existence proof one has to ensure that maximizers do in fact exist and that the 
shape constraint does not contribute to the Euler-Lagrange equation for the maximizer. 
With respect to the latter issue we introduce the improvement operator 

A®{AV) 


MV) := (1 - 5 ) 


\\AV{AV)\\, 


We are now able to describe the key arguments in the variational existence proof for 
solitary waves with profiles in C. 


Lemma 4 (three ingredients). 

(i) Since $ is strictly super-quadratic, each maximizing sequence for (If) is strongly 
compact. 

(ii) The cone C is invariant under the actions of both the convolution operator A and 
the superposition operator In particular, V E C implies R E C and %(V) E C. 

(Hi) We have V(Ts(V)) > V[V), where the equality sign holds if an only if V is a fixed 
point ofTs- 


Sketch of the proof. The main steps can be summarized as follows: 

(i) The assertion follows by a variant of the Concentration Compactness Principle. 

(ii) The invariance properties can be checked by straight forward calculations. 

(iii) Since $ is convex, we have 

V(Ts(V)) - V(V) > (d v V(Y), MV) - V) 

= a(y){%(v), MV)-V) 

= HV)\\MV)-V\\l, 

where we used ||75(E)|| 2 = ||E|| 2 = 1 — 5 and that cr(V) := ||A < h , (AI/)|| 2 /(l — 5) is 
well defined as long as V(V) > 0. 

The details can be found in [6]. □ 

Corollary 5 (variational existence proof). For any 0 < 5 < 1, there exists a solitary 
wave with Vs,Rs E C D BC°°(M) and ||V^|| = 1 — 5. Moreover, we have V{Vf) —* 00 as 
5 —»■ 0. 
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Proof. The existence of a solution Vs to (14) can be established by the Direct Method due 
to the compactness result from Lemma 4 and since V is strongly continuous. Moreover, 
any maximizer Vs satisfies 

V{V S ) > ViTsiVs )), ViVs) > V((l -S) X ). 

The first estimate implies V(Vs) = V(Ts(Vs)), so Vs satisfies the traveling wave equation 
(6) and is therefore smooth. We also compute 

V((l — 8)x) — 2 f $((1 — <5)(1 — x)) dx ——> +oo 

Jo 

and the proof is complete. □ 

The improvement operator Ts can also be used to compute solitary waves 
numerically. In fact, imposing homogeneous Dirichlet boundary conditions on a large 
but bounded and fine grid, the integral operator A can easily be discretized by Riemann 
sums. The resulting recursive scheme exhibits very good convergence properties; it has 
been applied to a wide range of potentials, see for instance [1, 6], and also been used to 
compute the numerical data displayed in Figures 1 and 2. 

A different variational framework has been introduced in [5] and later been applied 
to the high-energy limit in [3]. The key idea there is to minimize the kinetic energy 
term i||R||| subject to a prescribed value of p = V(V). The results from [3] imply 
for the potential (3) that solitary waves converge as p —> oo to the limit function 
X and satisfy Assumption 1 (though, strictly speaking, neither the unimodality nor 
the evenness of the profile functions R and V have been shown). In this context we 
emphasize again that uniqueness of Hamiltonian lattice waves is a notoriously difficult 
and an almost completely open problem. It is commonly believed that all variational 
and non-variational approaches provide - up to reparametrizations and for, say, convex 
potentials - the same family of solitary wave but there seems to be no proof so far. 

A non-variational existence proof for solitary waves with high energy has been given 
in [12] using a carefully designed fixed-point argument for the (negative) distance profile 
R in the space of exponentially decaying functions. In our notations, the smallness 
parameter is £ and the waves are shown to satisfy ||i? — Ax|| = 0 (e) in some suitably 
chosen norm. We therefore expect that the waves constructed in [12] also satisfy 
Assumption 1, although the justification of the unimodality might be an issue. 

2. Main result on the tip scaling of Rs 

Our first goal is to describe the asymptotic behavior of the distance profile Rs near 
x = 0 by showing that it converges as 5 — y 0 under an appropriately defined rescaling to 
some nontrivial limit function. In view of theorem 2 and the numerical simulations from 
figure 2 we expect that both the variable x as well as the shifted amplitude variable 
1 — Rs(5) must be scaled with certain powers of 5. A naive ansatz, however, does not 
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Figure 5. Cartoon of the tip scaling: The function Sg and its derivative for S > 0 
(gray, solid) and 5 = 0 (black, dashed). The dotted vertical lines enclose the symmetric 
interval Jg which has length l/yg ^ 1- The convergence Sg — » Sq as 5 —» 0 implies 
yg ~ £<s, see Theorem 7 and Corollary 8. 


work here because we lack a priori scaling relations between the small parameters S , 
and ps- Fc> r instance, if we would start with the rescaling 

R s (x) = 1 - S" n R 6 (5 12 x) , 

we could not eliminate as in the leading order equation. To overcome this problem 
we base our analysis on an implicit scaling, which magnifies the amplitude with Es but 
defines the rescaled space variable by 

x — fig x . 

In this way we obtain an explicit leading order equation that does not involve any 
unknown parameter and can hence be solved. Moreover, the corresponding solution 
finally allows us to identify the scaling relations between the different parameters; at 
the end it turns out that S, Es and ps are all proportional to each other, see Corollary 
13. 


2.1. Implicit rescaling of Rs 

In order to derive asymptotic formulas for Rs near x = 0, we define the rescaled distance 
profile 


Ss(x) := 


-Rs(O) - Rs(psx) 1 - £ 5 - Rsinsx) 


es 


es 


and obtain an even function which satisfies 

Ss(0) = $(0) = 0, S s (x) = Ss(-x ), 
see Figure 5 for an illustration. We also introduce the auxiliary functions 

hii) 

0,<x) :=eJ ,+1 4'(b 5 (-1+ W c)), 


(15) 

(16) 


as well as the intervals 


1 


( 17 ) 

(18) 
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Employing the identity (l)i as well as (4) we readily verify 

Sg(x) = — (Vg (p s x - |) - Vg {)igx + i)) 

= °8 £ ? +l ( v s {TsX - I) - Vg (figx + §)) 
and by (1) 2 we arrive at 

S'f(x) = 2 F s (x) - Gg(x) - Gg(-x ), (19) 

where we used that Rg(l + pgx) = Rs(~ 1 — Psx). The definition of Gs combined with 
the unimodality of Rs implies 

0 < G s {x) < G S (^j = e™ +1 &(R s (l)) < Ce™ +1 for all x e Jg (20) 
and from (17) we get 

T (eg + egSg(x) 


F s (x) = 


m +1 ’ 


( 21 ) 


^1 + S s (x) 

where the function T : [0, 1] —> M + with 


T(s) := $'(1 — s)s m+1 for 0 < s < 1, T(0) := lim T(s) =- 

s\o m + 1 

is smooth and positive. In particular, on the interval Jg we hnd 

2 1 


S%(x) « 2F s (x) 


™ +i h+s #( xr +i ’ 


and conclude that Sg satisfies the initial value problem (5) from the introduction up to 
small error terms. 

Lemma 6 (solution of the limit problem). The initial value problem (5) has a unique 
solution which is even, nonnegative, and convex. This solution S 0 grows linearly for 
x ±oo as it satisfies 


S' 0 (x) - psgn(x) 


< 


and 


xS' 0 (x) — S 0 (x) — K 
for all x eR with 

2 

p ■ = 


c 

(l + x) m 

(22) 

c 

(23) 

< 

— i ~\m—I 

(1 + X) 


r»oo 


yjm(m + 1) ’ 
and some constant C which depends only on m. 


k = x Sq(x) dx , rj / So(x) S'f (x) dir (24) 
Jo Jo 
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p, k, and q versus m 



Figure 6. Numerical values for the constants Ji, k, and rj from (24), which provide 
the leading and the next-to-leading order terms in the scaling relations between £$, ys, 
and ag, see Corollary (13). 


Proof. The planar and autonomous Hamiltonian ODE (5)i admits the conserved 
quantity 

2 

^pot(-t') ) E 'pot('l') ■ 


Etot(x) := 2 [S 0 (x, 


+ Er, 


\P 2 


1 + Sq(x) 


with value E to t(x) = .E tot ( 0 ) = \jl 2 for all x. A simple phase plane analysis reveals that 
S is even and that both S 0 and S' 0 are strictly increasing for x > 0, see Figure 4 for an 
illustration. In particular, we have 

S’ 0 (x) ^°° > \/2E tot (0) = JI > 0 

so the conservation law implies 

S’ 0 (x) - V / 2E'tot(0) 


2E'tot(0) — 2E pot (x) — s/2E tot (0) 
C 


< CE pot (x) < 


a . ~\m ? 

+ X) 

and hence (22). Moreover, the even function K 0 with K 0 (x) := xS' 0 (x) — 5^o(^) satisfies 

for all x > 0 , 


0 < Kq(x) = xSq(x) < 


C 


(l + x) m 

so K' 0 is integrable due to m > 1. The constant k = lim. 7,_>. O0 K Q (x) is therefore well- 
defined and (23) follows immediately from the estimate for K' 0 {x). Finally, rj is well- 
defined since the integrand is continuous and decays as x~ m for x —> oo. □ 

There seems to be no simple way to compute the constants k and rj as functions of 
m but numerical values are presented in Figure 6 . 


2.2. Asymptotic formulas for S ,5 

We now able to formulate and prove our main asymptotic result. 

Theorem 7 (asymptotics of Ss) • Any function S§ is strictly increasing and convex on 
Is- Moreover, the estimates 


sup 

x£Js 


Ss(x) - S 0 (x) 


< Ceff- 1 


(25) 
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sup 

S's(x) - S' 0 (x) 

< Ceff, sup 

Sf(x) - S"(x) 

xGJs 


x£Js 



< Ceff +l 


(26) 


hold for all 0 < 6 < 1 and a constant C independent of 6. 


Proof. Since S$ and S 0 are even functions it suffices to consider x E I, 5 . 
Dynamics of Ss- By (21) and due to T(0) — T(s) = Cs m+1 we have 


F s (x) - 


*( 0 ) 


m +1 


< Ce ™ +1 , 


1 + S s (x)^ 

and in view of (21) and (20) we conclude that S$ satisfies on the interval Is the ODE 


<§?(*) = 


vh(O) 


+ £™ +l h 5 {x) with \h s {x)\<C. 


m +1 


1 + S s (x)j 

Standard ODE arguments now imply 

S s {x) - S 0 {x) + S' s (x) - S' 0 (x) 


sup 

0 <o;<a;» 


< Ceff +1 


(27) 


(28) 


for any fixed ai* > 0, where C depends on 7*. 

Properties of Ss'- The monotonicity of both S§ and S' s on Is follows directly from 
the unimodality of Rs, Vs and the definition (15). In particular, Ss is convex on Is. We 
therefore have 

Sg(x) > S s (x *) + S s (x*) r ( x - xf ), 

and choosing ai* > 0 sufficiently close to 0 we find a constant c such that 

1 + Sg(x) > c(l + x) (29) 

holds for all x G Is, where we used that lim^ 0 S' s (x*) = S' 0 (x *) > 0 is implied by (28). 
Notice that (29) holds also for 5 = 0. 

Estimates for Ss- By (5) and (27) - and using both the Mean Value Theorem as 
well as (29) - we obtain 

S"(£) - So'(x) < C ]m +2 

(1 + x) 

Integration with respect to x yields 


Sg(x) - Sq(x) 


+ Ce? +1 


(30) 


S's(x) - Sq(x) 


< Ce^ +1 x + 


(i + y) 


c 

m +2 


S'Cz) - S' 0 (~Z) 


d5 d y 


since (16) ensures that 

Ss(0) = S o (0) = 0, S' s (0) = S' o (0)=0, 

and a direct computation reveals 


S' 5 (x) - S’ 0 (x) 


< Ce^ +1 x + 


(1 + 5 ) 


m +1 


S'sCz) ~ S' 0 (z) 


dz. 
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S's(x) ~ S' 0 (x) 


m+1 : 


< Ce'f 1 ‘xexp 


C 


(1 + 5 ) 


m+1 


dz < Ce™ +1 x , (31) 


and using x < 1 / (2//^) as well as the lower bound for ps from Lemma 3 we arrive at 
(26) i. Moreover, integrating (31) with respect to x gives 

5 s (x)-5 0 (x)| <Ce™ +1 x 2 

and hence (25). In combination with (30) we further get 


Sf(x) - S”(x) 
which in turn provides ( 26 ) 2 . 


< Ce™ +1 


x 


(1 + x) 


771+2 


+ 1 


□ 


A first consequence of Theorem 7 are leading order expressions for ps and hence 
for oy; below we improve this result by specifying the next-to-leading order corrections 
in Corollary 13. 

Corollary 8 (convergence of p$ and as). 

^5 <5—>o 

£5 


+ lb 




S —>0 _2 

-> IJ 2 


(32) 


Proof. By construction - see (15) - and Lemma 3 we have 
£ s Ss 


2 p 5 ) 1 55 R& 

and in view of Theorem 7 we conclude that 

1 \ <s-> 0 


5^-0 | 

-► I 


2e$ S 0 


‘dps J 


-> 1 . 


On the other hand, the estimates (22) and (23) evaluated at x = 1/(2 ps) imply 

(5-s-O 


2 ps Sq 


2 Vs) 




so the claim for ps follows immediately. The convergence result for as is a consequence 
of (4). □ 


3. Further asymptotic formulas 

In this section we exploit the asymptotic results on the tip scaling and derive 
approximation formulas for the foot and transition scaling, the tails of the profile 
functions, and for the scaling relations between 6 , £s, and as. In this way we obtain 
a complete set of asymptotic formulas which finally allows us to extract all relevant 
information on the limit 6 —> 0 from the function Sq only. We also sketch a possible 
application of these formulas, namely the study of the linearized eigenvalue problem. 
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1 
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/ 

1 

1 

X 
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j-U_ <1 > 

2w O UV) 


Figure 7. Cartoon of the tip and the foot scaling: The functions W$ and T$ for <5 > 0 
(gray, solid) and <5 = 0 (black, dashed). The limits Wq and T 0 can be computed from 
So, see (34) and (45). 

3.1. Transition scaling of Vs 

To describe the jump-like behavior of Vs near x = ±| we introduce the rescaled profiles 
Ws by 

Ws(x):=^V s (-± + ii S x) (33) 

£5 

and refer to Figure 7 for an illustration. The key observation for the asymptotics of Wg 
is the approximation 

W' s {x) « F s (x) « \~S'f(x ), 

and hence we are able to prove the convergence of Wg to 

W 0 (x) := i(^o(£) + /z) , (34) 

where Sq is the ODE solution from the tip scaling and defined in (5). 

Theorem 9 (convergence under the transition scaling). We have 

sup \W s (x) - W 0 (x) | < Cef , sup | W' s {x) - W'(x) \ < Ce™ +1 (35) 

xGJg xGJs 

for some constant C independent of 5. 

Proof. Uniform estimates for the derivative : By (33), the traveling wave equation (1), 
and the definition of fig in (4) we have 

w' s (x) = + mx) 

= Gg£™ +1 Vf (—| + flgx) 

= e^ +l ($'{Rs{nsx)) - &(Rg (-1 + figx))j , 
and from (17), (18) as well as (19) we infer that 

W' s ix) = F 5 (x) - G s (x) = i(^(£) - Gg(x) + G*(-x)) . 

The estimate (35)2 is now a direct consequence of (20) and Theorem 7. 
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Pointwise estimate at x — 1/(2/^): Using (33), the traveling wave equation (6), 


and the unimodality of Rs we obtain 


El 

hs 


W s 


£) = y ‘ (0) = 


■2 $'(i? 5 (x)) 


dx 


<?8 


.m +1 .*./ 


<^? +1 


f'(Ri(x)) dx , 


Thanks to (4), (17), and (19) we thus conclude 


W s 


M = 


1*8 


2 p s J as £™ +2 Jig 


2 F s (x)dx = / 2 F & {x)dx 


>h 


= / Sg(x) - G s (x) - G, 5 (-x) dx 

Jh 

and since f T Gs(Fx) dx = 0(e™) holds according to (20) and Corollary 8, we obtain 

W s \ 




2^5 ) 


^) +0(£ ”>- 


Lemma 6 combined with Theorem 7 and (32) finally provides 

H^) =p+o{£?)= H^) + o(£?) ’ 

so (35)i follows from (35) 2 . 


□ 


3.2. Tail estimates 

We complement our previous results by estimates for the tails of the profiles Vs and 
Rs- The derivation of those exploits the exponential decay with respect to x, which we 
establish by adapting an idea from [7]. 

Theorem 10 (tail estimates for Vs and Rs). The estimates 


sup Vs (x) 
W>i 


'M>i 


V s {x) dx < Cef 


and 


sup Rs (x) + 


hl>§ 


' W>§ 


R s (x) dx < Cef 


hold for some constant C independent of 5. 


Proof. Since Vs is unimodal and nonnegative by Assumption 1, we have 
0 < Vs (x) < Vs{— 1) for all x < —1. 

Moreover, the properties of the convolution operator A combined with (6) imply 
Rs(x) < Up + |)) K;f x\ < — - — for ail x < — 1 


(36) 


( 37 ) 
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and we infer that 

&(Vx(x + 1')') 

V s (x) < { 5K < Ce s m V 5 (x + 1) for all x < -2 (38) 

0’S 

where we used Corollary 8 , the properties of 4>', and that the unimodality of If 
guarantees 

2V S (-1) 2 < \\V 5 \\ 2 2 = (1 - 5f and hence V^(-l) < 1/^2 < 1 . 

By iteration of (38) we obtain 

V${x — n) < (Cs™) n Vs(x) for all x < — 1, n G N 
and conclude that Vs decays exponentially with rate 
A s >m |In o<51 (1 + o(l)) 


and satisfies 


-l 


V s (x) dx<CV s (-l). 


(39) 


Finally, (33) and Theorem 9 ensure that 


Ks(-l) = —w s , - — 
fJ-s V 2 ps 


£s_ 

As 


A + % 


2 As 


V 


\ 


0{e 


(40) 


= 0(e 


/ 


where the last identity stems from (22) and (32). The desired estimates for Vs are 
now direct consequences of (36), (39), and (40), and imply the claim for Rs thanks to 
(37) l □ 

Notice that the exponential decay rate in the proof of Theorem 10 is asymptotically 
optimal for small S. In fact, after linearization of 4)' in 0 we find the tail identity 

« A 2 V S , 


and the usual exponential ansatz predicts that the exact decay rate A^ is the positive 
solution to the transcendental equation 

sinh (A 5 / 2 ) _ —m /2 f 1 , f\ 

-- = y/ai = A£ S ^l + 0(e a )J. 

In particular, A s is large for small 5 and satisfies A a = m |lne a | + 0(ln |lne a |). 
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In this section we identify the scaling relations between the small quantities 


18 


(5, E$ , fJ.fi 

and start with two auxiliary results. 

Lemma 11 (scaling relation between p$ an d £5). The formula 
_ Jt £ s 
^ 1 £ s(f 1) 

holds for all 0 < 6 < 1. 


cyr 1 ) 


Proof. Our starting point is the identity 


1 - £<5 = Rs( 0) = 2 [ Vs(x) 
7-1/2 


dx 


- 1/(2 W ) 


— 2 ps 


V s [ -- + fix ) dx = 2e s 


"1/(2 Vs) 


Wsix) dx 


which follows from (6), Assumption 1, and (33). Theorem 9 now yields 

rVVus) 

Wq{x) dx + 0{ef) 


1 - Es = 2 Es 


£5 ( 


1^5 \ 2 

while (23) and (32) provide 


— — 7T + hsSo ( -— + 0{e 


1 \ 


2/1(5 / 


'<5 ) ) 


*(i) - (£) - ■* - + ° wo - K + ° wo' 

The combination of the latter two formulas yields 

l-£« = -fe-^«)+ 0(£j*), 

PS ' 1 

and rearranging terms we find the desired result. 

Lemma 12 (scaling relation between Es and 5). We have 


5 = 


1-J(l + £ s( K ~ 1)) (l — £ s (l + -jij j + 0{e™) , 


□ 


where the right hand side is real-valued and positive for all sufficiently small 8 > 0. 
Proof. From Assumption 1 and the tail estimates in Theorem 10 we derive 

(1 - 8 f = [ V 6 {xf dx = 2 f V 5 {xf dx + 0{e?) , 

7k 7-i 


Vs[xf dx = p 5 V 5 


i-i 


Js 


\ F 2 f 

-b psx ) dx = — / Ws(x) 2 dx . 

2 / p& J 

Js 


and (33) gives 
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Therefore 


and thanks to Theorem 9 - we get 


(1 - 5f = 2 ^ / W 0 {xf dx + 0(0 
Ps 


Js 


77 2 p 2 

V £ s 


+ 


2 V 5 VS Jl 


,§'(f) 2 dx + 0(O 


(41) 


where we used that S ' 0 is an odd function, and integration by parts yields 





thanks to the estimates and decay results from Lemma 6. In summary we find 

(l-i) 2 = ^-d(KJZ + i7)+0( £ J-) 

Vs VS 

and eliminating p$ by Lemma 11 yields via 

(1 — 5) = ^1 + £<5(ft — 1)^ — £^^1 + £^(ft — 1)^ ^ft + —^ + 0(£™) (42) 

the assertion. □ 


Lemma 11 and Lemma 12 provide explicit formulas for 
ps ~ £5 and 5 ~ £5 

in terms of £g but there is a slight mismatch between both results since the error bounds 
in the formula for ps are of higher order than those in the scaling law for 5. It is not 
clear, at least to the authors, whether this mismatch concerns the real error terms or 
just means that the bounds in Lemma 12 are less optimal than those in Lemma 11. 

Our main result concerning the scaling relations between the different parameters 
can now be formulated as follows. 


Corollary 13 (leading order scaling laws). The relation between p$ and £5 can be 
computed up to error terms of order 0(£™ +1 ) ? while the scaling law between 5 and 
£$ is determined up to order only. In particular, we have 


fti = HEi + II (1 — ft) E 2 t + o(eJ) , 

"s = PV + 2?Z 2 (1 - K) Ej m+l + o(ej m+1 ) 
for all m > 1, as well as 

2p — JIk — rj p 2 k 2 + rf + 2pKfj 2 


S = 


provided that m > 2 . 


2 p 


-£5 + 


’ v 2 


+ °{ £2 s) 


Proof. All assertions are provided by Lemma 11, Lemma 12, and formula (4). 


□ 
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3-4- Foot scaling of Rs 

We study now the asymptotic behavior of Rs near x — ±1. To this end we define 

Rs{ 1 + Msx) 


T s (x) := 


£5 


and find by direct calculations the identity 

T'\x) = F s (x) + Hs{x) - 2Gs{x) « \~S'f{x) 
because both Gg and 

H s (x) := &(Rs (-2 + fisx)) 
can be neglected on the interval Is- We further define 


T 0 (x) := !( S'o(x) + px + k 


and show that Ts converges as S —> 0 to T 0 , see Figure 7. 


(43) 

(44) 


(45) 


Lemma 14 (asymptotics of Rs at x = ±| and for Vs at x = 0). The terms R'g{— 
2 Rs(ts), and 14(0) are identical to leading order in 6. More precisely, we have 


R's{\) ~ + ^( £ T) — (1 + £ s( K ~ 1)) + 0( £ ™) j 


Ms 


and 


|2fl s (i) - = 0(ef) , |U(0) - = 0(e?) 

for all 0 < 5 < 1. 

Proof. From (15), Lemma 6, Theorem 7, and (32) we infer 

£s 




Ms 


2 Ms 


£s 


Ms 


2 Ms 


= -s' 0 [^)+out) = — + o(sT) 


Ms 


and similarly 


2 Ms 


Rs (!) = 1 - £s - £sSs ( ~— ) = 1 - £5 - £sS 0 ( 75 — ) +Q(m. 


2 Ms 


*?) 


= 1 -£s-£s \ 7 — - k + 0(e\ 


1 \ 


2 Ms \2 ms J 
= 1 + ej(R - 1) - ^ + 0(en • 

2 Ms 

Thanks to (33) and Theorem 9 we also hnd 




14(0) = 

1^8 


f~) + 0(e f) = ^£ + 0(e ?), 
2 Ms J Ms 


and the result follows from Lemma 11. 


□ 
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Theorem 15 (convergence under the foot scaling). The estimate 

sup (s 2 \T,(x) - f„(i)| + e I t' s (x) - T'(i)| + I Ti'(x) - T''(i)|) < Ce m+1 
xeJ s v J 

holds with some constant C independent of 5. In particular, we have 
R s {±l) = ±Ke s + 0{e?). 

Proof. The unimodality of R$, the monotonicity of < 3> / , and (20) imply 
0 < H s (x) < G s {x) < Ce ™ +1 

for all x £ J 5 . Combining this with (44) and Theorem 7 we arrive at the desired 
estimates for the second derivatives. We also notice that Lemma 6 along with Lemma 
14 imply 

T "(i) = ^“( 2 ^) +fr=-p+0(e7) 


= ^) +0 ( £ n = f;(- 


0 (e 


?), 


where the last identity stems from (43), and by similar arguments we justify 
1 


2/1,5 


- X 9 
2 S °\2ps 


1 
/' 

4(lH + hH* +0(£rl) 

- £ + Ae7A 

=+ o(C“‘) = f ‘(^) + . 

The assertions for the first and zeroth derivatives can thus be derived from the estimates 
for the second derivatives by integration with respect to x. □ 


1 Jt 1_ 

iw + 2 K 


3.5. Summary on the asymptotic analysis 
We finally combine all partial results as follows. 

Theorem 16 (global approximation in the high-energy limit). The formulas 


and 


£ 


Re(x) := < 


£ To 
0 


eS,[ ^ 

Te 


\x\ 


\x\ 


pe 


KM - J- 

pE 


0 


for 0 < |x| < \ 

for \ < |x| < | 
else 


(§ : M) 

for 0 < x 

V Pe ) 

else 
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with 


Pe ■ 


pe 


_-m— 2 -2 

' Pe 


1 + £ (k ~ 1 ) 

approximate the solitary waves from Assumption 1 in the sense of 


llife - 4,11, + \\v 6 - 4,11, + eTk, - & es \ = 0 (eT) = 0 ( 8 m ) 

for any q E [1, oo]. Here, S 0 solves the ODE initial value problem (5), the constants p, 
k are given in (24), and the functions W 0 , T 0 are defined in (34), (45)- 


Proof. Notice that 0(5) = 0(e, 5 ) is ensured by Lemma 12 and that it suffices to consider 
the case q — 00 because the functions R e and V £ are compactly supported. Theorems 7, 
9, and 15 - which concern the convergence under the different rescalings — as well as 
the tail estimates from Theorem 10 provide a variant of the desired estimates in which 
p £s is replaced by ps- Thanks to Lemma 11 we also have p$ ~ 5 as well as 

Ps = Pe s + 0(5 m+1 ) and hence — = + 0(5 m ~ 1 ). (46) 

pS p£ S 


Since S' 0 is bounded, the intermediate value theorem implies 


£,5 



£5 


Halloo M o{5 m - 1 ) 


0(5 m ), 


and by similar arguments we derive the corresponding estimate for T 0 . For the 
approximation of the velocity profile, the crucial estimate is 


Wn 


\x\ 


Ps 


bFn 


\x\ 


Pes 


ib«)6-M)o(V'- 1 ) = 0 (n. 


where £ denotes an intermediate value and where we used that the function x 1 —y xW((x) 
is bounded. Finally, the estimates for as — d £s follow from (4) and (46). □ 


For practical purposes it might be more convenient to regard e as the independent 
parameter and 5 as the derived quantity. In this case we can employ the following result, 
which is, however, weaker than Theorem 16 since the guaranteed error bounds are of 
lower order. 

Corollary 17 (variant of the global approximation result). We have 

II A, - fij.ll, + \\V, - UJI , + £"■-% - <T*,| = o( £ "->) 

for any q E [1, 00 ], where S £ := 1 — || 4||2 ~ £• 


Proof. Let e* be fixed, where the subscript * has been introduced for the sake of clarity 
only, and write <5* := S £t as well as p* := p Et for the quantities that can be computed 
directly and explicitly from e*. Lemma 12 provides 
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where £$, and ps, are defined by the exact wave data (Rs,, Vs,, cr s ,) and must not be 
confused with £* and /i*. On the other hand, a direct calculation - we just repeat all 
arguments between (41) and (42) with (e*, p,) instead of (e^, psf) ~ reveals 

+l/(2/n) 

(1 — S,) 2 = HKJI 2 = 2— [ W 0 (x) 2 dx 

P* J 

—1/(2m*) 

1 " e ’( 1 + l )) +0(er) ' 

Equating the right hand sides in both identities we then conclude 
£ s, = £* + 0(e™) 7 Vs, = p* + 0(e™) , 

where the last identity holds due to the ^-dependence of p * and Lemma 11, which 
provides an approximation of ps, in terms of £$,. Finally, exploiting the properties of 
So, Wo, and T 0 as in the proof of Theorem 16 we arrive at 

II 4. - Hoc + IIK, - KJU = o( k ‘~ £iJ ) = O {e ?- 1 ), 

and obtain analogous estimates for the g-norms due to the compactness of the supports. 
The assertion is now provided by Theorem 16. □ 

3.6. On the asymptotic eigenvalue problem 

Of particular interest in the analysis of solitary waves is the spectrum of the linearized 
equation. The problem consists of finding eigenpairs (A, U) G R x L 2 (M) such that 

XU = L 5 U , L S U := AQsAU , Q s (x) := , (47) 

where the function Qs becomes singular in the limit S —> 0, see Figure 8. Due to the 
shift symmetry of the traveling wave equation (6), there is always the solution 

A = l, U — Vf, 

and a natural question is whether this eigenspace is simple or not. In fact, simplicity 
would immediately imply some local uniqueness for solitary waves and is also an 
important ingredient for both linearized and orbital stability. 

Unfortunately, very little is known about the solution set of (47) due to the 
nonlocality of the operator A. In the small-energy limit of FPU-type chains, the 
corresponding problem has been solved in [4] by showing that the spectral properties of 
the analogue to Ls are governed by an asymptotic ODE problem which stems from the 
KdV equation and admits explicit solutions. The hope is that the asymptotic formulas 
derived in this paper provide spectral control in the high-energy limit. A detailed 
study of the singular perturbation problem (47) is beyond the scope of this paper but 
preliminary investigations indicate that the spectrum of Ls depends in the limit <5 —>■ 0 
- and at least for sufficiently large m - crucially on the coefficients c± 1 that are derived 
in the following result. 
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Figure 8 . Cartoon of the coefficient function Qs in the linear eigenproblem (47). 


Theorem 18 (weak*-expansion of Qs). For any sufficiently regular test function ip we 
have 

[ Q s (x)<p(x) dx = c_i/xjV(0) + c + i/^V'(0) + o(nf m{m ~ 2,3} 


with 


c-1 ■■= 


m+2 


dx, 


x 


C +1 • 2 


1 + Sq(x) 


1 + So(x) 
where the error terms depend on ip. 

Proof. Due to $"(r) = (1 — r)~ m 2 and the scaling relations (32) we find 

Qs(x)(p(x) dx = 0{af l ) = 0(p ™), 


m+2 


dx, 


J R\[-l/2,+1/2] 

and (15) along with (4) implies 


hs 


, Q s (x)<p(x)dx=^ I *» 

'[- 1 / 2 ,+ 1 / 2 ] a <5 Jj s 

1 f 


<f >"(R 5 (llsx))(p(p5x) dx 

p(hsx) 


^ (l + S 5 (x)) 


m+2 


dx. 


Theorem 7 as well as the linear growth of So - see Lemma 6 - ensure 


p(hsx) 


Js (1 + Ss(x)) 


m+2 


dx = 


p(psx) 


Js (1 + S 0 (x) 
<p(fj>sx) 


m+2 


dx + O^" 1 ) 


m+2 


1 + S 0 (x)J 

and by smoothness of <p and evenness of S 0 we can approximate 

f tp(nsx) 


dx + 0(/i™- 1 ) 


m+2 


dx = c_i</?(0) + c + i/i 2 +(0) + 0(t4) ■ 


1 + So(x) 

The claim now follows by combining all partial estimates from above. 


(48) 


□ 


For completeness we mention that the estimate (48) is not optimal for moderate 
values of m and might be improved for the prize of more technical effort. 
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